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Abstract 



> 

. We show how a number of NP-complete as well as NP-hard problems can be reduced 

to the Sturm-Liouville eigenvalue problem in the quantum setting with queries. We 
consider power queries which are derived from the propagator of a system evolving with 
a Hamiltonian obtained from the discretization of the Sturm-Liouville operator. We 
show that the number of power queries as well the number of qubits needed to solve 
the problems studied in this paper is a low degree polynomial. The implementation of 
power queries by a polynomial number of elementary quantum gates is an open issue. 
If this problem is solved positively for the power queries used for the Sturm-Liouville 
eigenvalue problem then a quantum computer would be a very powerful computation 
device allowing us to solve NP-complete problems in polynomial time. 



1 Introduction 

An important question in quantum computing is whether NP-complete problems can be 
solved in polynomial time, see E] and papers cited there. We address this question by 
studying the quantum setting with queries. We consider two types of queries: bit and power 
queries, see jE] for general information about queries and quantum computation. Here we 
only mention that bit queries are used in Grover's search algorithm [9 . They allow us to 
obtain the values of Boolean functions [3], and the approximate values of real functions 
|10j . Moreover, we know that bit queries cannot be used to solve NP-complete problems in 
polynomial time [3 . Power queries are used in the well-known phase estimation algorithm, 
see which plays a central role in Shor's factorization algorithm [21]. In a recent paper [20] 
we dealt with power queries in the study of the quantum complexity of the Sturm-Liouville 
eigenvalue problem. 

In this paper, we show how to reduce NP-complete problems to the Sturm-Liouville 
eigenvalue problem whose complexity in the classical and quantum settings has been studied 
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in [20]. Obviously, it would be enough to show this reduction for one NP-complete problem. 
We choose to present this reduction for several problems to show how the number of power 
queries and qubits depends on the particular NP-complete problem. In particular, that is 
why we consider satisfiability and the traveling salesman problem, as well as their NP-hard 
versions. The reductions presented in this paper can be summarized in the following diagram. 

SAT BOOL INT SLE 

TSP MIN BOOL INT SLE 

GRO BOOL INT SLE 

Here, SAT stands for the satisfiability problem, TSP for the traveling salesman problem, 
MIN for the minimization problem of choosing the smallest number out of N real numbers, 
GRO for Grover's problem, BOOL for the Boolean mean problem, INT for the integration 
problem, and finally SLE for the Sturm-Liouville eigenvalue problem. 

These reductions mean, in particular, that the satisfiability problem is reduced to the 
Boolean mean problem for a specific Boolean function which is reduced to the integration 
problem for a specific integrand, which in turn is reduced to the Sturm-Liouville eigen- 
value problem for a specific function, and finally the last problem is solved by the quantum 
algorithm using power queries. 

The Sturm-Liouville problem is defined in the next section. For the moment we mention 
that we want to approximate the smallest eigenvalue of a specific differential operator, and 
this smallest eigenvalue is given in a variational form as the minimum of specific integrals. 
We use a formula relating the Sturm-Liouville eigenvalue problem to a weighted integration 
problem, see [2D! ■ Many computational problems including the discrete problems mentioned 
above can be recasted as this weighted integration problem. Thus, we can solve them using 
the algorithms of [20] for solving the Sturm-Liouville eigenvalue problem. These algorithms 
use of order e -1 / 3 bit queries or log e _1 power queries and compute an ^-approximation of 
the smallest eigenvalue with probability |. The bounds on bit and power queries are sharp 
up to multiplicative constants, see 01201 • Hence, exponentially fewer power queries than 
bit queries are needed to solve the Sturm-Liouville eigenvalue problem. As we shall see, the 
same is true for the problems studied in this paper. 

In the quantum setting with bit queries, we do not obtain surprising results. The poly- 
nomial number of bit queries, e^ 1 ^ 3 , implies that the solution of NP-complete problems 
by modifications of the algorithm for the Sturm-Liouville eigenvalue problem will require 
exponentially many queries in terms of the NP problem size. 

The situation is quite different if we consider power queries. The logarithmic number of 
power queries, log e -1 , implies that NP-complete problems can be solved by modifications of 
the algorithm for the Sturm-Liouville eigenvalue problem and the number of power queries 
is polynomial in the problem size. 

More specifically, the satisfiability problem for Boolean functions with n variables can be 
solved with probability 1 — 5 using of order n log power queries and n qubits. Further- 
more, a truth assignment to a non-zero Boolean function with n variables can be computed 
with probability 1 — 5 using of order n 2 (log S^ 1 + log n) power queries and n qubits. 

The traveling salesman problem with m cities can be solved with probability 1 — 5 using 
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of order 



m log m(log 5 1 + log m + log <i max ) (log m + log d„ 



power queries and m log m qubits, where d max denotes the maximal distance between cities. 
Furthermore, an optimal route for the traveling salesman problem can be computed with 
probability 1 — 5 using of order 

m 2 log 2 m(log <5 -1 + log m) + mlog m(log 5' 1 + log m + log d max )(log m + log d max ) 

power queries and m log m qubits. 

Finally, Grover's problem for Boolean functions with n variables can be solved with 
probability 1 — 5 using of order n log <5 _1 power queries and n qubits. It is well known that 
Grover's problem requires of order 2 n / 2 bit queries. Therefore, it is evident that exponentially 
fewer power than bit queries are required for this problem. 

We stress that we only show how many power queries are needed to solve a particular 
problem. We use power queries which are of the form controlled- W Pj for a k x k unitary 
matrix W and some exponents pj. In our case, the matrix W is given by 

J¥ = exp(iiM g ) withi = v /z T. (1) 

Here, the matrix M q has a particularly simple form since it is a k x k real symmetric 
tridiagonal matrix, 



M q = (k + 1) 



2 -1 
-1 2 -1 



-1 2 -1 
-1 2 



+ 



vim) 



fc+i- 



with a function q : [0, 1] — > [0, 1] that is two times continuously differentiable and bounded 
by one up to the second derivative. This matrix corresponds to a classical approximation of 
the Sturm-Liouville operator. 

Contrary to the situation in Shor's algorithm where powers of a unitary operator can 
be implemented efficiently, the quantum implementation of the power queries for the k x k 
matrix W of the form Q by a number of known elementary quantum gates which is polylog 
in k is an open issue. If it turns out that the implementation cost of such power queries 
is disproportionally large compared to that of bit queries, then the positive results on the 
number of power queries are of only theoretical interest. If, on the other hand, power queries 
for W of the form ((IJ) and a function q satisfying conditions that we discuss later in this 
paper, can be implemented efficiently, i.e., at cost which is polylog in k, then we will have a 
very powerful computational device allowing us to solve NP-complete problems in polynomial 
time. 
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2 Sturm-Liouville and Integration 



We briefly recall the problem and some of the results from (2D]- We consider the following 
class of functions 



Q = \q : [0, 1] - [0, 1] 



qeC 2 ([0,l]) and max ||g (i) ||oo <1 

where C 2 ([0,1]) stands for the class of twice continuously differentiable functions, and 
Halloo = max^o,!] \q(x)\. 

We studied the approximate computation of the Sturm-Liouville smallest eigenvalue X(q) 
which is defined in the variational form by 

. Jo[(u'(x)) 2 + q(x)u 2 (x)] dx 

X{q) = nun -y . 9 . , , , (2) 

o^ueH^ J Q u 2 (x) dx 

where Hi is the Sobolev space of absolutely continuous functions for which u' G -^([0, 1]) 
and u(0) = u(l) = 0, see [21 [22]- Combining results from from UHl 123] we have the 
formula that relates the Sturm-Liouville smallest eigenvalue problem to integration, 

f (q(x) - |) sin 2 (vrx) dx = \ (X(q) - vr 2 - §) + O (\\q - fQ . (3) 
Jo 

We analyzed the quantum setting with bit and power queries in j2D] . For bit queries, we 
showed that \(q) can be computed with error r\ and probability | using 9(r] _1//3 ) bit queries, 
and this bound is sharp modulo a multiplicative constant. It is easy to check that from this 
result follows that NP-complete problems of size n can be solved with an exponential in n 
number of bit queries. Therefore, from now on, we restrict ourselves to the quantum setting 
with power queries for matrices W of the form (0). 

In j2Hl we presented a quantum algorithm (p based on phase estimation applied to the 
discretized matrix of the Sturm-Liouville problem. The initial state was an approximate 
eigenvector, as proposed by Abrams and Lloyd 0, and computed by the algorithm of the 
Jaksch and Papageorgiou [14J. The algorithm (p computes X(q,rj) such that 

|A(g) — \(q, rj)\ < r] with probability | Vg G Q 

using of order log77 _1 power queries, log 2 ?? -1 additional quantum operations, 0(1) function 
values of q and classical operations, as well as log?] -1 qubits. 

It is well known that we can increase the probability of success to, say, 1 — 5 by repeating 
the algorithm of order log 5 _1 times and then taking the median as the final approximation. 
The algorithm <p with repetitions computes X(q,r),S) which is an //-approximation of X(q) 
with probability 1—5, i.e., 

\X(q) - X(q,i],5)\ < r] with probability 1-5 Vg G Q. (4) 

The resulting algorithm with repetitions cf) uses of order 
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• \ogS~ 1 log ry - 1 power queries, 

• \og5~ 1 log 2 ?] -1 additional quantum operations, 

• 0(1) function values of q and classical operations, and 

• log??" 1 qubits. 

In the next sections we show how to modify the algorithm <p to solve a number of con- 
tinuous and discrete problems. In what follows, we restrict ourselves and mention only the 
number of power queries and the number of qubits of these modifications because they are 
the most important characteristic of the cost of a quantum algorithm. The rest of the cost 
characteristics can be easily derived from the corresponding components of the cost of the 
algorithm with repetitions. 

We start with integration. Knowing how to approximate X(q) we can approximate the 
integral in (J3j) modulo the second term which is of order \\q — ^H;^,. We provide the details 
in the next section. 



3 Integration 

Consider the (weighted) integration problem 

/(/) := / f(x) sin 2 (7nc) dx 
Jo 

for functions / from the class 

Fa, = {/eC 2 ([0,l]): max ll/^IU <M}. 

Here M is a positive number. We want to compute an ^-approximation of /(/) with proba- 
bility 1 — 5 on a quantum computer with power queries. Since |/(/)| < M we assume that 
e < M since otherwise is an ^-approximation, and the problem is trivial. Without loss of 
generality, we also assume that e < 1. 

Observe that / G implies that the function 

Qf,d x ) = I + c f( x ) Vx e [ML 

belongs to Q, defined in the previous section, for c G (0, (2M) -1 ]. In this case, the formula 
(J3J) states 

Define 

e e 2 

c = ~rTn~i t and ri = ce = — — — T . 

M 2 log E" 1 ' M 2 log e- 1 
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Let \(qf tC ,r),5) be an ^-approximation of A(g/ )C ) with probability 1 — 5 computed by the 
algorithm with repetitions of the previous section. Knowing \(qf, c , V> $) we compute on a 
classical computer 

Then 

\I(f) - A Int (f,e,5)\ < —v + 0(cM 2 ) = ±e(l + o(l)) with probability 1-5. 

Hence, for small e, A lnt (f, e,5) is an ^-approximation of /(/) with probability 1 — 5. In 
this way we can solve the integration problem and we summarize this result in the following 
theorem. 

Theorem 3.1. We compute an e- approximation with probability 1 — 5 for the integration 
problem for the class ¥ M by the quantum algorithm A lnt using of order 

• log^" 1 (log M + log e" 1 ) power queries and 

• log M + log e~ l qubits. 



4 Preliminaries 

We now present some preliminaries that will be used as technical tools to translate the inte- 
gration problem of the previous section to the NP-complete and NP-hard problems discussed 
in this paper. 

Take a function h G C 2 ([0, 1]) with h®(0) = = for i = 0, 1, 2, and for which the 

integral h(x) dx is positive. Examples of such functions include h(x) = (x(l — x)) a with 
a > 2 or h(x) = x 3 (l — x) 3 g(x) for a positive g G C 2 ([0, 1]). We extend the domain of h by 
defining H(x) = h(x) for x G [0, 1] and H(x) = otherwise. Due to the boundary conditions 
imposed on h, we have H G C 2 (M). 

For a positive (large) integer N, we subdivide the interval [|, |] by introducing the points 

x i = \ + Ijf for J = 0,l,...,iV. 
For j — 0, 1, . . . , N — 1, we define the functions 

h ^ = W~2 H ( 2N ( X ~ x i)) for x e [°» !]■ 

Observe that hj vanishes outside the open interval (xj,Xj+i), and ||/ij||oo = ||^||oo/(4iV 2 ), 
Halloo = \\h'\\J/(2N), and Halloo = Halloo. Hence, if we set 

M: = — (ii?.^.ii A "ii~). 
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which is equal to ||/i"||oq for large N, then 

hj GF M for j = 0,1,...,JV-1. 

Observe finally that 



where 



/ hj(x) dx = / hj(x) dx = — — Int(/i) for j = 0, 1, . . . , N — 1, (5) 
Jo Jxj 8N 

Int(/i) := / h{x)dx. 
Jo 



5 Boolean mean 

Consider the class B n of all Boolean functions of n variables mapping {0, l} n into {0, 1}. We 
can equivalently assume that the domain of such Boolean functions is {0, 1, . . . , N — 1} with 
N = T. Hence, BgB„ means that 

B:{0,l,...,iV-l} -> {0,1}. 

We want to approximate the mean 



3=0 



by a quantum algorithm with power queries. 

We now show how this problem can be reduced to the integration problem of Section |H1 
Using the notation of Section HJ we define the function / B : [0, 1] — > E by 



hj(x) B(j) (2 sin 2 (7rx)) 1 if x G [xj,x j+1 ] with j = 0, 1, . . . , N - 1, 
ifxG [0 i]U[|,l]. 



Observe that due to the fact that hj vanishes up to the second derivatives at Xj and Xj + i, 
and the fact that 2 sin 2 (7ra) > 1 for x G [j, |], we conclude that /b G C 2 ([0, 1]). Furthermore 

II/bHoo < ||/ i ||oo/(4iV 2 ), 

Halloo < ||/i'||oo/(2iV)(l + 0(iV- 1 )),a^ 

ii/i'iu < ii^iua+o^- 1 ))- 

Hence, if we set M = ||/i"||oo(l + 0(A^ -1 )) then / B G F M . Observe that, due to ©, the 
integration problem for takes now the form 

T ,r, v^i™, x T J+1 , / x , Int(/i) 1 %=$ m , , Int(/i) n 

3=0 ' A ' j=0 



Let T) = lnt(h)e/ (16iV 2 ). We now use the quantum algorithm A lnt (f B , rj, S) from Sectional 
which, for small e or large N, computes an ^-approximation of I(fm) with probability 1 — 5. 
Knowing A(f®,r},5) we compute on a classical computer 

1 fi/V 2 



Then 



1 fi N 2 1 fi N 2 

muom-^ oo1 (B,M)| = -^—\I(f M )-A^(f B , v ,5)\ <-4ttt^ 



Int(fc) 1 - M(/l y 

and this holds with probability 1 — 5. We summarize this result in the following theorem. 

Theorem 5.1. We compute an e- approximation with probability 1 — 5 for the Boolean mean 
problem for the class M n by the quantum algorithm A^°° l using of order 

• log5 _1 (log e^ 1 + n) power queries and 

• log e~ l + n qubits. 



It is known that the amplitude amplification algorithm of Brassard, H0yer, Mosca and 
Tapp [Hj computes an e-approximation of Sat(B) with probability 8/tt 2 = 0.81... using 
of order min(iV, e~ l ) bit queries. Furthermore, this number of queries is order-minimal as 
proven by Nayak and Wu ^Hj- We stress that the basic part of the quantum algorithm A^° o1 
is the phase estimation algorithm which uses power queries. For e > iV -1 this algorithm 
uses of order log e" 1 power queries to compute an ^-approximation with probability 8/7T 2 . 
Hence, the nunber of queries has an exponential improvement in its the dependence on e^ 1 . 

The amplitude amplification algorithm of 6J has been used as a basic tool for solving 
many continuous problems such as real mean, multivariate integration, path integration and 
multivariate approximation in the quantum setting. These problems have been defined over 
many classical spaces such as L p , Sobolev and Korobov spaces. For a number of these 
continuous problems, the bit query complexity and the quantum speedups over the worst 
case and randomized settings have been established based on the optimality of the quantum 
summation algorithm, see jH3 HH [TJJ HH1 UH1 EH Ej The use of power queries yields an 
exponential improvement in the number of queries. This can be achieved simply by using 
the exponentially better power query bound for the Boolean mean. 



6 Satisfiability 

The satisfiability problem, SAT for short, is a well known NP-complete problem in the 
Turing machine model of computation |Sj. This means that all NP-complete problems can 
be reduced to SAT in polynomial time, and if the conjecture P^NP is true then there are 
no algorithms solving SAT in polynomial time with respect to the length of the Boolean 
function expressed in conjunctive normal form. SAT can be stated as a decision or as a 
computational problem and we deal with both of them in this section. 



As in the previous section, consider the class B n of all Boolean functions of n variables 
with the domain {0, 1, . . . , N — 1}, where N = 2™. The two variants of SAT problem are 
defined as 

• SAT X : for IB G B n given in the conjunctive normal form, verify if there exists an index 
j such that B(j) = 1. 

• SAT 2 : for a non-zero B G B n given in the conjunctive normal form, compute an index 
j such that B(j) = 1. 

We now show that either problem can be solved with probability 1 — 5 by using a number 
of power queries which is polynomial in n and log 

We begin with SATx and use the notation of Section 03 Observe that for any B G B n , 
the mean Sat(B) is a multiple of A -1 , i.e., Sat(B) = k/N for some G {0,1,..., N}. If we 
have a real number x such that |Sjv(B) — x\ < ^ then \k — N x\ < \ which implies that 



Obviously, k > iff there exists an index j for which B(j) = 1. 

For r] < we conclude that from an 77-approximation of Sjv(B), with probability 
1 — 5, we can compute the exact value of S(B) with probability 1 — 5. We know that, for 
small//, j4 Boo1 (B, rj, 5) computes an //-approximation of Sat(B) with probability 1 — 5. Letting 
rj = 1/(3N) we compute A® ool (B, 1/(3N), 5), which is an l/(3A^)-approximation of S N (M). 
From this we can compute the exact value of Sat(B) with probability 1 — 5. This means that 



solves the satisfiability problem. We stress that to compute A^ AT (B, 5) we run the quantum 
algorithm A Bool (B, l/(3iV), 5) and the rest is computed on a classical computer. Since we 
know how many power queries and qubits are used by v4 B ° o1 we obtain the following theorem. 

Theorem 6.1. The satisfiability problem SAT\ for the class B n is solved with probability 
1 — 5 by the quantum algorithm ^ ATl which uses of order 

• n log<5 -1 power queries and 



We turn to the SAT 2 problem. That is, for a non-zero B from B n we want to compute 
an index j G {0, 1, . . . , N - 1} for which B(j) = 1. 

We will use bisection on the domain of B. Every bisection step will shrink the cardinality 
of the domain by 2. Using the quantum algorithm ^ ATl with k — n, n — 1, . . . , 0, we will 
know whether an index of the true assignment belongs to the the decreased domain. In this 
way, after n steps we identify an index j for which B(j) = 1. Since A^. ATl is a probabilistic 



k = [Nx + |J and S 7V (B) 




N 




• n qubits. 
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algorithm and we use it n times, we need the success probability of this algorithm to be 
1 — Si, where 

(l-<fi) n = 1-5. 

For small 5, we obviously have Si = S/n(l + o(l)). 
More precisely, let 

D k = {0, 1, . . . , 2 k - 1} for k = 0, 1, . . . , n. 
We set j n — 0, and perform the following steps for the Boolean function B from B n . 
For k — n — 1, n — 2, . . . , 1, do: 

• define f k : D k ^ {0, 1} by f k (j) = B(j + j k+1 ) for j E D k , 

• run the quantum algorithm A k iATl (f k ,Si) and compute on a classical computer 



3k = 



Jk+ i + 2 k ifAf Tl (/ fe A) = NO, 
Jfc+1 HA s k AT ^f k ,Si) = YES. 



Finally we set 

A^(M,S) = Jo . 

We claim that the algorithm A n 2 solves the satisfiability problem SAT2, i.e., for the index jo 
we have B(j ) = 1- 

Indeed, first note that j k < 2 k + 2 k+1 + ■ ■ ■ + 2 n ~ l < 2 n - 1 = iV - 1 and j + j k+1 < 
2 k — 1 + jk+i < N — 1 for j G -Dfe. Therefore, the Boolean functions are well defined. 

For the first step, k — n — 1, we have /„_! = B on the first half, D n _i, of the domain D n . 
We check whether f n _i is zero. This holds with probability 1 - Si iff A^ifn-i, Si) = NO. 
If fn-i = then B is non-zero on the complement of D n _i and an index j for which B(j) = 1 
is at least equal to 2 n ~ 1 . That is why we define j n -i = 2 n ~ 1 in this case. If, however, / n _i 
is non-zero over D n _i then we are looking for an index j with B(j) = 1 in the set D n -i, and 
we set j n _i = j n = 0. In this way, after the first step we can restrict the search of an index 
to the domain of cardinality 2 n ~ l . For the second step it is enough to work with the domain 
-D„_2 and use the proper shift j n -i in the definition of the Boolean function f n -2- After n 
steps we identify a proper index j = j for which B(j) = 1. In fact, it is easy to see that we 
will find the smallest index j for which B(j) = 1 

Since the quantum algorithm A^ ATl (f k , Si) works with probability 1 — Si and we repeat 
n times the algorithm, the probability of success is at least (1 — Si) n . By the definition of 
Si, this is equal to 1 — S. This proves the following theorem. 

Theorem 6.2. The satisfiability problem SAT2 for the class B n is solved with probability 
1 — S by the quantum algorithm Af^ 2 which uses of order 

• n 2 (log5 _1 + log n) power queries and 

• n qubits. 
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7 Grover's problem 



Grover's problem can be defined as the satisfiability problem SAT2 for the class M n of Boolean 
functions of n variables for which we know a priori that there exists exactly one index j = j& 
for which B(j) = 1. Obviously, the algorithm Af/^ 2 solves Grover's problem although the a 
priori knowledge about the uniqueness of the index j B is not used. We now show that using 
this a priori knowledge it is possible to find a more efficient quantum algorithm than ^ AT2 . 
Define the weighted Boolean mean 



W N (M) 



1 N ^ 



N 



for 



3=0 



Clearly, Wjy(B) = jn/N. Hence, it is enough to compute the exact value of Wn(M) and then 
j'b = NWn(M). This can be achieved by switching to the integration problem, as we did in 
Section for the function : [0, 1] — > R defined by 



j hj(x) B(j) (2 N sin 2 (7rx)) if ] with j = 0, 1, 





N-h 



ifx€ [0,J]U[|,1]. 



As in Section we conclude that 

In 

1(9*) 

Hence, 



2N 



e ¥ M for M 

x m+ 1 



\h"\ 



0(N~ 1 )), and 



hj (x) dx 



Int(/i) 
16iV 3 



W 



N 



1QN 4 
Int(/i) 



In Section EJ we defined the algorithm A lnt such that, for small e, A lnt (f, e, 5) is an 
er-approximation of /(/) with probability 1 — S. 

Let e = Int(/i)/(48iV 4 ). Define the quantum algorithm 

16iV4 Ah*, n 1 



Grover / 



Then 



16iV 4 
Int(/i) 



.4 



5) 



Int(/i) 



16iV 4 
lnt(h) 



I(g M )-A™(g M ,e,5) < 



1QN 4 
Int(/i) 



1 

3' 



and this holds with probability 1 — 5. Hence, 

j B = Af ovcr (B,5) with probability 1 - 5. 
This and Theorem 13.11 yield the following theorem. 

Theorem 7.1. Grover's problem for the class B n is solved with probability 1 — 5 by the 
quantum algorithm v4^ rover which uses of order 

• n log power queries and 

• n qubits. 
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8 Minimization 



In this section we consider a real number minimization problem which we will use to solve 
the traveling salesman problem in the quantum setting with power queries. 
For positive N = 2 n and M, define the set 

X n ,M = {x = [x , xi, . . . , xat_x] : Xj G R and \xj\ < M for j = 0, 1, . . . , N — 1 }. 

Let 

Mm(x) = min x~. 

i=0,l,...,JV-l J 

The minimization problem is defined as: 

• MINi : compute A(x) which is an e-approximation of Min(a;) with probability 1 — 5, 
i.e., |Min(x) — A{x)\ < s holds with probability 1 — 5 for all x G ^ n ,M- 

• MIN 2 : compute an index j = j(x) for which |Min(x) —Xj\ < e with probability 1 — 5 
for all x G X n M . 

Clearly, using a classical computer we must use each Xj at least once and that is why 
the worst case and randomized complexities are proportional to N, i.e., they are exponential 
in n. We now show how to solve this problem in the quantum setting using a number of 
power queries which is polynomial in n, log M, log 5" 1 and log e^ 1 . 

We begin with the minimization problem MINi. For a real number y, define the Boolean 
function f y : {0, 1, . . . , JV - 1} -»• {0, 1} by 



fy(j) 

Then let 



1 ifxj < y, 
if Xj > y. 



1 N ^ 

3=0 

Note that 

S N (f y ) > iff y > Min(x). 

Clearly, the condition Sjv(/ y ) > is equivalent to the SATi problem for the Boolean func- 
tion f y and can be solved by the quantum algorithm A^ ATl of Section |3 

Initially, we know that Min(x) G [—M,M]. That is why we set ao = —M, bo = M and 
2/o = 0, and use the bisection algorithm for the interval [— M, M\ with k* steps, 

■ M 

We also choose 5\ such that (1 — 5\) k * = 1 — 5. For small 5, we have 6\ = 6/k*(l + o(l)). 
More precisely, we perform the following steps. 
For k = 1, 2, . . . , k* do: 
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• run the quantum algorithm A SATl (f yk _ 1 , Si), 



• compute on a classical computer 




Vk-i 
h-i 



Vk-i 



if^Ok-^i) 



if^Uk-xA) 



YES, 
NO, 



YES, 
NO, 



• and y k = \(a k + b k ). 



Finally set 



After fc* steps we have the interval [a k *,b k *] of length 2M/2 k * < 2e, and Min(a;) G 
Therefore |Min(a;) — < e is an approximation of Min(x). Note that this 
algorithm works with probability (1 — <5i) fc * = 1 — 5. Knowing the requirements of the 
algorithm A^ ATl we obtain the following theorem. 

Theorem 8.1. We compute an e- approximation with probability 1 — 5 for the minimization 
problem MIN\ for the class X ni M by the quantum algorithm A^ 1 ^ 1 using of order 

• n (log M + log e -1 ) (log<5 -1 + log (log M + log power queries and 



We now turn to the minimization problem MIN2. It is easy to see that this problem can 
be solved by combining the quantum algorithms developed so far. To explain the main idea 
of the quantum algorithm for the minimization problem MIN2 we ignore for a moment the 
fact that all the quantum algorithms of the previous sections work probabilistically. Knowing 
V — Vk* by the A^ ini algorithm, we apply the A^°° l algorithm of Section |3] and compute the 
exact value of Sjv(/y)- If Sjv(/y) 7^ then f y is a non-zero Boolean function. If S(f y ) = then 
y < Min(x) and since y is an ^-approximation to Min(x), we have Min(x) —s<y< Min(x). 
In this case we have Min(x) < y + e and < y + e — Min(x) < e. Hence, y + e is also an 
^-approximation to Min(a;) and S(f y+£ ) 7^ 0. Thus we modify y :— y + e if S(f y ) = 0. Then 
S(fy) 7^ in either case, and y is still an ^-approximation to Min(x). Knowing that f y is 
a non-zero Boolean function, we now run the A^ AT2 algorithm and compute an index j for 
which f y (j) = 1, or equivalently, Xj < y. Hence, Min(x) < Xj < y < Min(x) + e and Xj is 
an ^-approximation that solves the minimization problem MIN 2 . We now formalize the idea 
of this algorithm. Let (1 — 5\) 3 = 1 — 5. Hence, for small 5 we have 5\ = 35(1 + o(l)). We 
perform the following steps: 

• run the quantum algorithm A^^(x, e, 5i) to obtain y, 

• run the quantum algorithm A® ool (f y , 1/(3 AT), 81) to obtain z, 

• if z = then set y := y + e, 



• n qubits. 
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• run the quantum algorithm A^ AT2 (f y , Si) to obtain j. 
Finally set, 

Then the index j solves the minimization problem MIN 2 and this algorithm works with 
probability (1 — 5i) 3 = 1 — 5. Counting the number of power queries and qubits of all parts 
of the algorithm we obtain the theorem. 

Theorem 8.2. We compute an e -approximation with probability 1 — 5 for the minimization 
problem MIN2 for the class \ n ,M by the quantum algorithm A^ 1 ^ 2 using of order 

• n 2 (log S~ l + log n) + n (log M + log e' 1 ) (log<5 -1 + log (log M + log 

power queries and 

• n qubits 



9 Traveling salesman 

The traveling salesman problem, TSP for short, is a well known NP-complete problem that 
deals with the shortest tour between m cities, see |H]. Let D = [d(j, k)]™ k=l denote the m x m 
matrix with d(j, k) being the distance between the city j and the city k. We assume that 
d(j, k) are positive integers for j ^ k and d(J,j) = 0. The TSP can be studied as a decision 
or as a computational problem and we consider the three variants: 

• TSPi : for a given integer B verify if there is a permutation ir = [71* (1), 7r(2), . . . , n(m)} 
of indices {1,2,..., m} for which 

m 

d(ir) := £d(7r0V0' + l)) < B 

with 7r(m + 1) = 7r(l). 

• TSP2 : compute 

Min(L>) = mind(7r). 

7T 

• TSP3 : compute a permutation 7i~* such that 

d(vr*) = Min(D). 

Observe that if we can solve TSP2 then it is just enough to check whether Min(Z^) < B. 
Similarly, if we can solve TSP3 then it is enough to compute d(n*) to solve TSP2. 

We now show that TSP is a special case of the minimization problem studied in the 
previous section. Indeed, consider the set P m of all m\ possible permutations, and let g : 
{0, 1, . . . , m\ — 1} — > P m be an injective mapping. We take 

n := [log ml] = mlog m (1 + o(l)). 
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Let iV = 2". For j = ml, ml + 1, . . . , N — 1 we extend the function g by setting g(j) 
g{ml - 1). Then g : {0, 1, . . . , JV - 1} -»• P m and 

g(j) — 7Tj = [7r i (l),7r i (2),...,7r i (m)] for j G {0, 1, . . . , JV - 1}. 

Defining 

m 

Xj = d(TTj(k),irj(k + 1)) with ^-(m+l) = TTj(l), 



fc=i 

we see that 



min Xa = Mm(D). 

j=0,l,...,N-l J 



Note that Xj can be computed using a number of bits which is polynomial in m and the 
maximal distance between cities, 

ci max = max d(i,j). 

i,j=l,2,...,m 

Furthermore, Xj G [0, M] for j G {0, 1, . . . , N — 1} with 

M > md max . 

To run the algorithms A^ 1 ^ and A n ^ for the solutions of the minimization problem we 
need to know an upper bound on m<i max . This can be achieved as follows. For k = 0, 1, ... , 
we run the quantum algorithm A^ ATl (l — / 2 fc ^i) °f Section El to check, with probability 
1 — Si, whether the Boolean function 1 — j 2 k defined in Section |H] is zero. If so, then / 2 fe = 1 
meaning that for all Xj < 2 k , and we can take M = 2 k . Hence, we perform as many steps 
as necessary until the first occurrence of A^ ATl (l — /2 fc ; £l) — NO. This will happen after at 
most p := [log 2 m + log 2 <J m ax] steps. Then we define M = 2 P and M G [mrf ma , 2mrf max ]. 
We choose Si such that (1 — 5i) p = (1 — S) 1 ^ 2 to conclude that we obtain M with probability 
(1 - S) 1 ' 2 . For small S, we have 8 X = S/(2p)(l + o(l)). 

Knowing M, we can solve TSP by solving the minimization problem of Section |HJ Since 
all Xj are integers, it is enough to compute an ^-approximation y of Min(x) with e < | to 
conclude that [y + |J is equal to Min(x). The number y can be computed with probability 
(1 — S) 1 / 2 = 1 — S 2 by the quantum algorithm A^ 1 ^ 1 (x, e, 5 2 ) with e, say, |. In the same way, 
we can find a shortest path represented by a permutation ix* by first computing the index j 
by the algorithm A^ 1 ^ of Section |H] for which |Min(a;) — Xj\ < |. Again, since Xj is an 
integer, Xj = Min(x) and 7r* = g(j) is a needed permutation. Hence the following algorithms 

. x_ if L<^M^) + U < 5 ' 




I NO if otherwise 

Al SP %D,S) = i^ (j i a) = 

solve the traveling salesman problems TSPj for j = 1, 2, 3, respectively, with probability 
1 — 5 2 = (1 — 5) 1 / 2 . We also know the proper M with probability (1 — S) 1 ^ 2 . Therefore, the 
whole algorithm works with probability 1 — 5. This proves the following theorem. 
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Theorem 9.1. The traveling salesman problems TSPj are solved with probability 1 — 5 by 
the quantum algorithms im Pj using of order 

• for j = I, 

— mlog to (log 5~ l + log to + log d max ) (log to + log d max ) power queries and 

— m log m qubits, 

• for j = 2,3, 

to 2 log 2 to (log 5' 1 + log to) 

+ TO log TO (log (T 1 + log TO + log rf max ) (log TO + log rf max ) 

power queries and 

— mlog to qubits. 
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